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Statistical Behavior of Rain Attenuation* 

By S. H. LIN 

(Manuscript received October 25, 1972) 

Thirty-one sets of experimental data on the statistics of microwave rain 
attenuation at frequencies above 10 GHz, in the U.S.A., England, Japan, 
Italy, and Canada, indicate that: (i) the distribution of rain attenuation a, 
in dB, is approximately lognormal with a standard deviation a a of logio a 
ranging from O.46 to 0.7 1 for earth-space paths, and from 0.33 to 0.86 
for terrestrial paths; (ii) the distribution of the rain fade duration r is also 
approximately lognormal with a standard deviation a T of logio t ranging 
from 0.44 t° 0. 76 for both earth-space paths and terrestrial paths. We 
propose a theory to explain this general behavior. A theoretical upper 
bound for the fade duration distribution in the tail region is also given. 

The findings in this paper simplify the determination of rain attenua- 
tion statistics needed for the design of earth-satellite radio links and ter- 
restrial radio links. 

I. INTRODUCTION 

The statistics of rain attenuation are important for the design of 
both terrestrial and earth-satellite radio links using frequencies above 
10 GHz. Many experiments 1-33 have been performed to obtain data 
on rain attenuation for different frequencies, path lengths, and geo- 
graphical locations. This paper presents the behavior of rain attenua- 
tion statistics found in our study to be common to the available ex- 
perimental data. 1 

II. DEFINITIONS 

Let 

V(t) be the time- varying amplitude of the received signal 
voltage normalized to its nonfaded level, 



* An excerpt of this paper has been presented at the 1972 IEEE International 
Conference on Communications at Philadelphia, and included in the Proceedings of 
the Conference. 

' These data are actual measured microwave rain attenuations and not the ap- 
proximate attenuations calculated from rain rate data. 
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a (t) = — 20 logio V(t) be the time- varying rain attenuation 
in dB, 
P(a ^ A) be the expected fraction of time that a(t) exceeds any 
specified value A, 
P be the expected fraction of time that rain falls at the 
location of the radio link, 
P c (a ^ A) be the expected fraction of raining time that a ex- 
ceeds A, 
a m be the median value of a during the raining time, i.e., 
Pc(a ^ a m ) = 0.5, 
t(A) be the duration of rain attenuation fades with a ex- 
ceeding any specified threshold A in dB, 
P[t(A) ^ b] be the probability that the fade duration t(A) ex- 
ceeds any specified duration b, and 
f ( A ) be the average duration. 

Notice that 

P(a ^ A) = Po-Pcia ^ A). (1) 

Thus, P may also be called the probability of rainfall and P c (a ^ A) 
the conditional distribution of a under the condition that the rain is 
falling. 

III. SUMMARY OF RESULTS 

3.1 A ttenuation Distribution 

(i) The available experimental data on both earth-space paths 
and terrestrial paths in the U.S.A., England, Japan, Italy, and 
Canada show consistently that the conditional distribution, 
P c (a ^ A), is approximately lognormal within the attenua- 
tion range, 1 dB ^ a ^ 50 dB, of practical interest.* 
(w) We propose a theory to explain the lognormal behavior of 
attenuation a(t). In essence, the value of the rain attenuation, 
a (t) , at any time instant can be multiplicatively affected by a 
large number of random time-varying parameters of the 
environment such as the present states and the past histories 
of the weather conditions at various locations all over the 
world. The large number of random multiplicative components 
and the central limit theorem lead to the lognormal distribu- 
tion of a (t) . 



* This means the distribution of signal amplitude 7(0 is approximately 
log-lognormal. 
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{Hi) These results indicate that three parameters, P , a a (the 
standard deviation), and a m (the median), are sufficient to 
determine the rain attenuation distribution P{a ^ .A). These 
parameters depend on geographic locations. 

(iv) a m increases almost linearly with the path length because the 
median rain rate usually is small and is almost uniform over the 
entire path. 
(v) The dependence of a m on frequency follows the theoretical 
prediction of Sezter 20 for small rain rate because the median 
rain rate is usually small. The rain rate data in Refs. 13 and 
34 to 37 indicate that the median rain rate ranges from 0.5 to 
10 mm/h, depending on location. 

(vi) a a decreases slightly as the path length increases because of 

the averaging effect of the propagation volume. 2425 
(vii) <r a decreases slightly as frequency increases. 

3 . 2 Fade Duration Distribution 

(i) The experimental data indicate that the fade duration dis- 
tributions, P\_t(A) ^ b~\, are also approximately lognormal. 
The physical reason for the lognormal duration is the same 
as that for the lognormal attenuation. 
(m) The probability of occurrence of long fade durations has an 
absolute upper bound given by 

pCfczW-^)' (2)* 



2 



for any X ^ 1. For example, P(r/f ^ 10) ^ 0.0161 means 
that no more than 1.61 percent of the total number of fades 
will have durations longer than 10 • f. 

IV. EXPERIMENTAL DATA 

4. 1 Rain A ttenuation Distribution 

Thirty-one sets of experimental data on rain attenuation distribu- 
tion, each with a time base of six months or longer, are summarized 
in Tables I, II, and III for earth-space paths, long terrestrial paths 
(> 10 km), and short terrestrial paths (^10 km), respectively. 1 " 

* The constant M is denned in eq. (13). 

1 In the literature, there are many other sets of experimental data with a time 
base less than six months. Those data are not included because the short-term dis- 
tributions of rain attenuation are fairly random. 
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(1) 30 GHz, NEW JERSEY 

(2) 16 GHz. NEW JERSEY 

(3) 19 GHz, ENGLAND 




RAIN ATTENUATION A IN dB 

Fig. 1 — Lognormal distributions of rain attenuation on earth-space paths. 



100 



The original data are all given in terms of the unconditional dis- 
tribution P(a ^ A). Based upon P (either available or estimated), 
we convert these unconditional distributions into the conditional dis- 
tribution P c (a ^ A) by eq. (1). This conversion is done to exclude the 
dry periods in which the rain attenuation is identically zero.* 

When plotted on a lognormal coordinate system, these conditional 
distributions, P c (a ^ A), are all approximately straight lines within 
the attenuation range, 1 dB ^ a ^ 50 dB, of practical interest. 
Figures 1 to 3 show nine examples. The equation describing the 
lognormal distribution is 



d / ^ a \ i f r logi ° A ~ m °i 

P c (a ^ A) = \ erfc 

L V2c7 J 



(3) 



where erfc (~) denotes the complementary error function; <r a is the 
standard deviation of logm a during the raining time ; and 



Ha = logioa, 



(4) 



is the mean value of logioa during the raining time. The estimated 
values of <r a and a m are given in Tables I, II, and III. 

However, the accuracy of these estimated values of a a and a m is 



* Absorption by the clear atmosphere is not considered in this paper. 
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1)11 GHz, 43.5 Km, ALABAMA 

(2) 12.62 GHz, 80 Km, JAPAN 

(3) 15 GHz, 15.78 Km, CANADA 

(4) 18 GHz, 24 Km, ENGLAND 




100 



RAIN ATTENUATION A IN dB 
Fig. 2 — Lognormal distributions of rain attenuation on long terrestrial paths. 

limited by two problems : 

(i) The time base of the experiment may not be sufficiently long 

to yield stable statistics, and 
(m) some of the published experimental data do not provide the 

probability P„ of rain during the experiment.* 

Let 

m CtogioafflH-M.. (5) 



If a is lognormally distributed, then (3 will be normally distributed with 
zero mean and unity standard deviation. This allows us to pool all 
the available data of on the same graph paper for comparison. Figures 
4, 5, and 6 show the pooled data of the earth-space paths, long ter- 
restrial paths, and short terrestrial paths, respectively. It is seen that 
these experimental results of /3 are indeed normally distributed with 
zero mean and unity standard deviation. 

* The ? values at Crawford Hill, New Jersey, are provided by D. C. Hogg and 
R. A. Desmond from their rain gauge records. The P„ values in Alabama, England, 
Japan, and Canada are estimated from the information in Refs. 6, 8, 11, and 13. 
The P„ values in Italy and Massachusetts are assumed values using some judgment. 
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(1) 18 GHz, 5 Km, GEORGIA 

(2) 18 GHz, 10 Km, ENGLAND 




100 



RAIN ATTENUATION A IN dB 

Fig. 3 — Lognormal distributions of rain attenuation on short terrestrial paths. 

However, Fig. 6 shows that the data of short paths have significant 
deviations from the normal distribution in the tail region (/3 ^ 2.8). 
Furthermore, most of the deviations in the tail region are downward 
from the straight line approximation. The reason for this nonsym- 
metric deviation is discussed in Appendix B. 



POOLED DATA OF SIX EXPERIMENTS 
ON EARTH - SPACE PATHS 
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Fig. 4 — Normal distribution of 8 of earth-space paths. 
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POOLED DATA OF FOURTEEN EXPERIMENTS ON LONG PATHS 




Fig. 5 — Normal distribution of /3 of long terrestrial paths. 



90 r- 



POOLED DATA OF ELEVEN EXPERIMENTS ON SHORT PATHS 




Fig. 6 — Normal distribution of of short terrestrial paths. 
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1 - 



(1) 2.6 Km, CRAWFORD HILL, N.J. 

(2) EARTH - SPACE, CRAWFORD HILL, N.J. 

(3) 20 Km, ITALY 

(4) 24 Km, ENGLAND 



(4) 




FREQUENCY IN GHz 

Fig. 7 — Effect of frequency on median attenuation a m . 



4. 2 Dependence of <j a and a m on Path Length and Frequency 

The experimental results, discussed in Section 4.1, show that the 
three parameters, P , a a , and a m are sufficient to determine the dis- 
tribution P(a ^ A). Therefore, it is important to study the depend- 
ence of (7 a and a m on path length and frequency. 

The effects of path length and frequency on a a and a m are shown in 
Figs. 7 to 10. These experimental results indicate the effects of path 



1.6i— 



1.0 - 



0.5 



(1) 2.6 Km, CRAWFORD HILL, N.J. 

(2) EARTH - SPACE, CRAWFORD HILL, N.J. 

(3) 20 Km. ITALY 

(4) 24 Km, ENGLAND 




20 



40 



FREQUENCY IN GHz 

Fig. 8 — Effect of frequency on standard deviation <r„. 



568 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1973 

Si- 
ll) 18.5 GHz, CRAWFORD HILL, N.J. 

(2) 17 GHz, MOBILE, ALABAMA 

(3) 11 GHz, ENGLAND 

6h 



? 4 - 




30 

PATH LENGTH IN Km 

Fig. 9 — Effect of path length on median attenuation a m . 

length and frequency on <r a and a m as stated in terms (iv) to (vii) of 
Section 3.1.* 

The theoretical calculation 20 indicates that the increase of rain at- 
tenuation with frequency / is slightly faster than the square law in the 
range : 10 GHz ^ / ^ 60 GHz, and 1 mm/h ^ rain rate ^ 15 mm/h. 
This is the basis of the dashed curves in Fig. 7. The slopes of straight 
lines in Fig. 9 are proportional to median rain rate and extinction 
coefficient. Curves (1) and (2) in Fig. 9 indicate that the median rain 
rate in Alabama is much larger than that in New Jersey. 

4.3 Fade Duration Distribution 

The available nine sets of experimental data on the histogram of the 
durations t(A) of rain attenuation fades are summarized in Table IV. f 
We convert these histograms into the cumulative distribution 
P[_t(A) ^ b~]. On a lognormal coordinate system, these fade dura- 
tion distributions are all approximately straight lines. Figure 1 1 shows 
two examples. 

* In Fig. 8, the a in England seems to increase slightly with frequency in contrast 
to those in the U.S.A. and Italy. A possible reason for this inconsistency is that the 
time bases for the three sets of data for 11, 18, and 36 GHz measured in England 
are not concurrent. 

+ Some of the available data in the literature are not included because of short 
time base. 
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1.6 r- 



1.0 



0.5 



(1) 18.5 GHz, CRAWFORD HILL, N.J. 

(2) 17 GHz, MOBILE, ALABAMA 

(3) 11 GHz, ENGLAND 



10 



20 



30 



40 



50 



60 



PATH LENGTH IN Km 
Fig. 10 — Effect of path length on standard deviation a a . 

The equation describing the lognormal distribution of t(A) is 

"lOglQ b — fir' 



P[t(A) ^ 6] = | erfc 



V2 



<7 r 



(6) 



where n T and o> are the mean and the standard deviation respectively 
of logio t(A). The estimated values of <x T and the average fade dura- 
tion f(A) are given in Table IV. Again, the unstable statistics, caused 
by insufficient time bases, limit the accuracy of these estimated values 
of <t t and f. 
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(1) 30 GHz, 21 dB FADE THRESHOLD, N.J. 

(2) 19 GHz, 5dB FADE THRESHOLD, ENGLAND 




FADE DURATION b IN MINUTES 

Fig. 11 — Lognormal distributions of fade duration. 
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V. THEORY 

5.1 Rain Attenuation Distribution 

Existing theory 19-28 and experimental data 1-19,29-33 indicate that 
rain attenuation is a complicated function of many parameters of the 
propagation medium : the total number of rain drops in the path, the 
drop size distribution, the fine grain spatial characteristics of the rain 
density along the path, wind velocity, the presence of up or down 
drafts, raindrop shape, raindrop cant angles, the storm cell shapes and 
sizes, raindrop temperature, etc. In other words, the attenuation a(t) 
is a function of many random time-varying parameters of the medium. 
Furthermore, through the coupling of the atmosphere, the above- 
mentioned parameters of the propagation medium depend on the 
present states and the past histories of the weather conditions at 
many near or faraway locations and altitudes. 

We will assume that the rain attenuation can be affected by a large 
number of random time-varying multiplicative components : 

a(t) = S l (t)-S 2 (t)-S 3 (t)---S n (t). (7) 

Each of [5<(0}'-i represents the random modification factor due to 
an environmental parameter. 

Taking logarithms on both sides of eq. (7) yields 

log a = log Si + log S 2 + ■ • • + log S„ (8) 

which shows that log a is a summation of a large number of random 
variables. Then by the central limit theorem, 38 the distribution of 
log a approaches a normal distribution for large n if there is no domi- 
nant component. Therefore, the distribution of a is approximately 
lognormal.* t 

The basis for the multiplicative formulation (7) is that the environ- 
mental parameters affect the rain attenuation a(t) in a proportional 
fashion (i.e., in terms of percentage) rather than an additive fashion. 
For example, at 30 GHz frequency, the theoretical calculation shows 
that, when the rain temperature decreases from 20°C to 5°C, the rain 
attenuation increases by approximately 4 percent. This means the 

* It is interesting to note that, in Fig. 12 of Ref. 10, the lognormal distribution 
appears to be a reasonable fit to the probability distribution of rainfall rate. The 
physical reason for this behavior probably is similar to that for rain attenuation a 
discussed in this section. 

* A basic characteristic of a lognormal random variable y is that its value can vary 
in the entire semi-infinite range : g y < w . For the rain attenuation problem, 
the attenuation a(dB) possesses this basic characteristic. On the other hand, the 
normalized signal amplitude V(t) is strictly confined to the finite range : ^ V(t) ^ 1, 
which rules out the possibility of a lognormal distribution for V(t). 
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variation of rain attenuation, due to 15°C variation of rain tempera- 
ture, can be 4 dB, 0.4 dB, or 0.04 dB if the attenuation at 5°C is 100 
dB, 10 dB, or 1 dB, respectively. Similar arguments apply to the 
effects of other parameters on the attenuation. A more general inter- 
pretation of formulation (7) is discussed in Appendix A. 

We emphasize that some of the components {<S<(0}<-1 may be 
extremely slowly varying functions, which may take several months 
or even several years in order to show their effects. For example, the 
rainfall intensity-duration-frequency data of the Weather Bureau 32 ' 39 
show that in New Jersey the return period 40 for a rain rate exceeding 
150 mm/h, which continues for a 5-minute duration, is about 5 years. 
Therefore, if the time base of a rain attenuation experiment is less 
than 5 years, the chance of missing these rare and extreme events is 
very high. The justification for considering such a long-term distribu- 
tion of a (t) is that the microwave radio systems, which are designed 
based on these statistics, contains many repeaters, each sampling its 
own rain universe and contributing to total path outage. 

5.2 Fade Duration Distribution 

We also assume that the duration t(A) of a rain attenuation fade 
with a ^ A is affected by a large number of random time varying 
multiplicative components 

t(A) = X 1 -X i -X,---X n . (9) 

Each of \Xi\iZi represents the random modification factor of an 
environmental parameter. Therefore, the long-term distribution of 
t(A) is also approximately lognormal. 

5.3 Upper Bound for Fade Duration Distribution 

In radio system design, one is concerned with the occurrence prob- 
ability of an unusually long continuous outage. It is desirable to have a 
"quick estimate" of the fade duration distribution, especially in the 
tail region of long duration. In our experimental and theoretical 
study 41-43 of the lognormal distributions of durations of rain and 
multipath fading, W. T. Barnett 41 has found an upper bound on the 
fade duration distribution as discussed in the following. 

Let* 

x = -■ (10) 

f 

* The idea of normalizing fade duration r to the average duration f originates 
from the work 48 of S. 0. Rice on fade duration distributions. 



RAIN ATTENUATION 



573 



The lognormal distribution of x can be written as 

"log™ X — n~ 



*«-*-p£^| 



(id 



where n and a are the mean and the standard deviation, respectively, 
of log™ x. For lognormally distributed x, it is easily shown that 

where x is the mean value of x, and 

M = log™ e a* 0.434. (13) 

The definition (10) implies that 



x = - = 1. 

f 



Equations (12) and (14) show that 



M = 



Substituting (15) into (11) yields 



P(x ^ X) 



= 5 erfc 



lcr 2 
2M' 

logic Z + a*/2M' 



(14) 



(15) 



(16) 



Therefore, the lognormal distribution of x is completely determined 
by only one parameter a. 



Pm«x (-=->10) = 1.61% 
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10 1.5 2.0 

STANDARD DEVIATION a 

Fig. 12 — Effect of standard deviation a on the probability of long duration of fade. 
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By differentiating (16) with respect to <r, it is easily shown that 
dP 
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• at <t = <2M log w X 



(17) 
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Pig. 13— Maximum probability for various normalized fade durations. 
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which implies that for any A r ^ 1 the probability P(x ^ I) as a 
function of a has a maximum value : 



Pm.*(x ^ X) = \ erfc 



at 



log™ XI 
M _ 



* = V2M log 10 X. 



(18) 



(19) 



For example, Fig. 12 shows P(x ^ 10) as a function of a. 

Figure 13 shows the maximum probability Pmax(x ^ X) as a func- 
tion of X as given by eq. (IS). By the use of "C-discriminant equa- 




NORMALIZED FADE DURATION X (= y ) 

Fig. 14 — The envelope of the family of lognormal duration distributions. 
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tion" in the theory of ordinary differential equation, 45 it can be shown 
that eq. (18) is, in fact, the envelope of the family of lognormal dis- 
tributions (16) with a as the family parameter. In other words, the 
family (16) are all tangent to (18) as shown in Fig. 14. The shaded 
area in Fig. 13 is a forbidden region where the lognormal duration dis- 
tribution will never penetrate. Therefore, the upper bound for the 
fade duration distribution is given by eq. (2). 

VI. CONCLUSION 

Both the experimental data and theory indicate that both the rain 
attenuation distribution and the fade duration distribution are 
lognormal. The detailed results have already been given in Section III 
(Summary of Results). 
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APPENDIX A 

Generalized Interpretation of Formulation ( 7) 

From a more general viewpoint, the effect of an environmental 
parameter Z t (t) on a(t) may be more complicated than the simple 
linear proportional relation 

a(t) cc Zi(t), i= 1,2, ••-,«. (20) 

For example, a(t) may depend not only on the present value, but also 
on the past history of the environmental parameter Z«(')i tnen tne 
relation between a(t) and Z,-(0 becomes 



a{t) cc / H^t, nZiit'W 

i = 1, 2, • • • , n 

where Hi(t, t') is the impulse response of a(t) if the input Zt(t) is an 
impulse 8(t - t') applied at t'. In order to allow for the more general 
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and complicated relation between a(t) and Z,(£), we shall use the 
mathematical operator notation: 

a(t) oc GlZi(t)l i = 1, 2, 3, ...,n (22) 

where G, is a transformation from Z,(0 into a(Q. Then the formulation 
becomes 

a (0 = <7i (Zi) ■ G 2 (Z 2 ) ■ ■ • (? n (Z„) . (23) 

Let 

5,(0 = G,[Z,(0] » = 1, 2, 3, • -, n. (24) 

Substituting (24) into (23) leads to (7). 

The formulation can be further generalized to include the cases 
where the effects of various environmental parameters on a are not 
completely separable. Let 

7 = logioa, (25) 

ah = logio Si = logio Gi(Zi) 

% - 1, 2, 3, •••, n - 1, and (26) 

R = Iogxo 5 n . (27) 

Substituting (25), (26), and (27) into (8) gives 

7(*i, h, ■ ■ ■, fn-l) = 7 + Cl(fc " 10 + C 2 (*2 - |») + • • • 

+ C n _ 1 (^ n _ 1 - |„_0 + (B - R) (28) 
where 

7 = Cih + C 2 | 2 + • • • + C»_i|„_, + ft. (29) 

Equation (28) can be interpreted as the first-order Taylor series ex- 
pansion of y(|i, £ 2 , ••■, £ n -i) with a remainder term,* ft — ft. The 
main point in eq. (28) is that the effects of the (n — 1) random vari- 
ables { t-i\\Zl~ x on 7 do not have to be completely separable because 
of the remainder term R — ft. If R — ft does not dominate the sum in 
the right-hand side of eq. (28), then, by central limit theorem, the 
distribution of y (i.e., logio a) is approximately normal even if the 
effects of various environmental parameters are not completely 
separable. (This argument originates from Refs. 46 and 47.) Therefore, 
the formulation (7) includes very general and complicated relation- 
ships between the environmental parameters and rain attenuation a (t) . 



* If higher-order derivatives of y(^, £ 2 , • • -, £„_i) exist, then R — R~ represents 
the sum of all the higher-order terms; otherwise, R — R~ represents the difference 
between y(|i, {2, ■ ■ •, £ n -i) and its first-order Taylor series expansion. 
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APPENDIX B 

Deviations of Experimental Data From Lognormal Distribution 

Aside from the experimental error, there are two major factors which 
contribute to the deviations of experimental data from the lognormal 
distribution. 

B.l Effect of Finite Number of Components 

Davenport and Root 38 have indicated that when the number, n, of 
components is finite, the normal distribution may well give a poor 
approximation to the tails of the distribution of the sum (8) even 
though the limiting form of the sum distribution is, in fact, normal. 
Therefore, in practice, we believe that the deviation of the experimental 
data from the lognormal distribution may increase toward the tails. 

B.2 Effect of Time Base 

Since some of the components in eq. (7) are extremely slowly 
varying, reducing the time base will reduce the number of contribut- 
ing components since the influences of slow components are approxi- 
mately constant in a short experiment. Therefore, we expect the 
deviation of the experimental data from the lognormal distribution to 
increase as the time base decreases. 

The minimum required time base for the convergence of the experi- 
mental data to the lognormal distribution increases as 

(i) path length decreases, or 
(ii) operating frequency decreases, or 
(Hi) attenuation range of interest increases. 

Furthermore, when the time base of a rain attenuation experiment 
is not long enough, the deviations of the short-term distribution 
Ps(a ^ A) from the long-term distribution P(a ^ A) are usually non- 
symmetric in the deep fade region, i.e., P s (a ^ A) is more likely to 
be less than P(a ^ A). This nonsymmetry is caused by: 

(i) The distribution of fade duration t(A) is lognormal, which is 
nonsymmetric with respect to the average fade duration 
f(A). Typically, about 70 percent of fade durations are 
shorter than the average duration. 
(ii) The probability distribution of the number N(A, T) of ob- 
served deep fades, exceeding the margin A dB in the period T, 
is somewhat similar to a Poisson distribution which is also 
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nonsymmetric with respect to the average number N(A, T) of 
deep fades in a period T. 

Since, when the time base is too short, both N(A, T) and t(A) have 
a higher chance of being less than their average values N(A, T) and 
f(A), respectively, then the short-term distribution P,{a ^ A) tends 
to deviate downward from the long-term distribution in the tail 
region.* 

B.3 Effect of Time Base on Standard Deviation a a 

Since it takes a long time base to include appreciable effects of 
slow components, we expect that <r a increases slightly as the time base 
increases, and reaches an asymptotic value only after the time base 
is long enough to include the effects of all the possible slow components. 
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